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Abstract 

We compute the 0(a t a s ) two-loop corrections to the neutral Higgs boson masses in the 
Minimal Supersymmetric extension of the Standard Model. An appropriate use of the effective 
potential allows us to obtain simple analytical formulae, valid for arbitrary values of and of 
the mass parameters in the stop sector. We elucidate some subtleties of the effective potential 
calculation, and find full agreement with the numerical output of the existing diagrammatic 
calculation. We discuss in detail the limit of heavy gluino. 
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1 Introduction 



There is a crucial prediction of the Minimal Super symmetric extension of the Standard Model, 
or MSSM (for reviews and references, see e.g. ref. Q), subject to decisive tests at present and 
future colliders. It is the existence of a light CP-even neutral Higgs boson, h, accompanied by 
other states (H, A, H^), whose masses are strongly correlated but can vary over a wide range 
near the weak scale. An important step in the understanding of the MSSM Higgs sector was the 
realization that the classical bound < mz, and, more generally, the classical relations among 
the gauge and Higgs boson masses, are violated by large radiative corrections, dominated by top 
and stop loops |||, ^ |||. After that, extensive efforts have been devoted to progressive refinements 
of the theoretical predictions for the Higgs boson masses and couplings, as functions of the 
relevant MSSM parameters. These activities have been performed in several directions, with 
special emphasis on the prediction for m^: inclusion of stop mixing effects ||, q|; resummation of 
large logarithms using appropriate one- and two-loop renormalization group equations (RGE) 
P, [?|, p| S|; complete one-loop diagrammatic calculations including momentum-dependent effects 
flO| , [TT[| ; calculations of the most important two-loop contributions [12, 1^, 14, 15]. Other studies 
have been oriented towards a meaningful combination of the above results [15, 16 1, and towards 



the implementation of the latter in computer codes JlT], 18], to be used in turn for experimental 



analyses [19]. 

In the present paper we address once more the computation of the neutral Higgs boson masses, 
whose present state-of-the-art can be summarized as follows. There is a diagrammatic two-loop 
computation, including (D(ata s ) effects, performed for arbitrary uia and arbitrary values of the 
parameters in the stop mass matrix, in the zero-momentum limit |13|] . While, for small stop 
mixing and universal soft stop masses, sufficiently simple and accurate analytical formulae have 
been obtained [14], in the general case the complete formulae are rather lenghty, which may be 
a problem for their practical implementation in computer codes. The results of ref. fl3| can be 
improved by including the logarithmic 0(af ) corrections, as extracted by solving perturbatively 
the appropriate RGE There is also a computation [15] of both 0(ata s ) and O(a^) two- 

loop corrections to m^, based on the effective potential approach. This computation, however, 
is applicable only for m^ 3> mz- Moreover, the full results of for nih are available only in 
numerical form, and accurate and simple analytical formulae were provided under the additional 
assumptions of small stop mixing and universal soft stop masses. 

In view of the situation described above, there is still room for a number of useful improve- 
ments that could be achieved without excessive effort. As a first step, one should aim at simple 
analytical formulae for the two-loop corrected mass matrix of the neutral CP-even Higgs sector, 
still in the zero-momentum limit and including only 0{atta s ) corrections, but for arbitrary values 
of rriA and of the parameters of the stop mass matrix. One could then proceed with the inclusion 
of the O(ol) corrections, and of the corrections coming from the momentum-dependent part 
of the two-loop Higgs propagators, into the above framework. Finally, one could address the 
resummation of the large logarithms of (m^/m^), the ratio of the two stop mass eigenvalues, 
by means of suitable RGE, defined in an appropriate effective theory. This has been done |§ in 
the case of small stop mixing, but is considerably more complicated in the case where a large 
splitting between m; and is induced by a large mixing term in the stop mass matrix. 

In this paper we accomplish the first step of the above program, leaving the remaining 
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steps for future work. The paper is organized as follows. After this introduction, section 2 
recalls the general features of the calculation of the MSSM neutral Higgs boson masses in the 
effective potential approach. Section 3 describes the main features of our two-loop calculation 
of the 0{ottcts) contributions, and presents its results, in a form that allows to assign the input 
parameters either in the DR scheme or in some on-shell scheme. In the concluding section we 
compare our results with the existing literature, and discuss in detail the heavy-gluino limit. 
This limit requires some care, especially when the input parameters are assigned in the DR 
scheme, as often done in models that predict the soft supersymmetry-breaking masses. 

Technical details are confined to three appendices. Appendix A gives the analytical expres- 
sions for the two-loop contributions to the neutral Higgs mass matrices that are controlled by 
the gluino mass. Appendix B gives the explicit formulae that are needed for the transition from 
the DR scheme to our implementation of the on-shell scheme. Appendix C gives the relation 
between m\ and m\, which may be useful for discussing models that predict the values of the soft 
supersymmetry-breaking masses, and in particular of mj, at some cut-off scale for the MSSM. 



2 Higgs masses in the effective potential approach 



The MSSM neutral Higgs boson masses can be identified with the zeros of the corresponding 
two-point functions, which depend on the external momentum and have in general the form of 
a matrix. In the limit of vanishing external momentum, these masses can be formally obtained 
by the following method: compute the effective potential V e s, retaining its complete dependence 
on the neutral Higgs fields, H® and H®; minimize V e g to find the vacuum expectation values 
(Hf) = vi/y/2 and (H®) = V2/V2; expand V e s around its minumum up to quadratic fluctuations, 
and diagonalize the resulting mass matrix. Of course, at each step we can carry out explicitly 
only the calculations that contribute to the final results, at the desired level of approximation. 
Before moving to our specific two-loop computation, we give now some formulae that illustrate 
the formal results of the general procedure, and are valid at every order in perturbation theory. 

Putting all the other fields to zero, and keeping only the dependence on the neutral Higgs 
fields, the tree-level Higgs potential of the MSSM reads: 



V 



where: m\ 
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Hi 



Hi 



+ m 2 2 



ml 
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m, 



+ m| (fljflj + h.c.) + (\H^ - \H^Y , (1) 



H 2 + H ' A 4 * s Higgs mass term in the superpotential; m Hi , 
fn 2 H2 and m\ are soft supersymmetry-breaking masses; g and g' are the SU (2) l and U (l)y gauge 
couplings, respectively. It is not restrictive to choose m\ real and negative, so that v\ and V2 are 
real and positive, and the neutral Higgs fields can be decomposed into their vacuum expectation 
values plus their CP-even and CP-odd fluctuations as follows: 



Hi = 



v 1 + S 1 + iP 1 



v 2 + S 2 + iP2 



(2) 



V2 ' z ~ V2 

In the effective potential approach, the mass matrices for the neutral CP-odd and CP-even 
Higgs bosons can be approximated, at every order in perturbation theory, by: 
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where V e g = Vq + V is the loop-corrected Higgs potential in the DR scheme, and (Si) = (Pi) = 
{S2) = {P2) = 0. Using the explicit expression of the tree-level potential, eq. (|T|), v\ and V2 are 
determined by minimizing the effective potential: 
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Combining eqs. ©-(H), the CP-odd and CP-even Higgs mass matrices become (i, j = 1, 2): 
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Combining further eqs. (^) and (^), we can write the CP-even Higgs mass matrix as follows: 



(Ml) ={-!)»* (M 2 P ). 
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where the CP-odd mass matrix can be expressed, in terms of the loop-corrected CP-odd Higgs 
mass rriA and of tan/3 = v^/vx, as: 

sin 2 (3 sin (3 cos j3 
, sin (3 cos cos 2 j3 
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m 2 A . 



(9) 



3 0{a t oi s ) two- loop corrections to the neutral Higgs boson masses 

The formulae derived in the previous section have general validity, and were employed long ago for 
the one-loop computation [g, H |B|, j| . We will now follow the same strategy for the calculation 
of the 0{a.tots) two-loop corrections to the entries of the neutral CP-even Higgs boson mass 
matrix. The relevant Feynman diagrams involve top, stop, gluons and gluinos on the internal 
lines, and are shown in figure |l]. Since V e g must be considered in a generic Higgs background, it 
is important to elucidate the dependence of the propagators and vertices on the Higgs fields. 

In a generic Higgs background, the MSSM Lagrangian contains the following bilinear terms 
in the top fields: 

=(t[t£)('t - x : ) (it), do. 



-x 



where t' L and t' R are four-component fermions of definite chirality, and the field-dependent mixing 
term X is: 

X = h t H% = \X\e iv , (0<v9<2vr). (11) 

It is not restrictive to assume that ht is real and positive, then = htV2/V2 and (93) = 0. 

Analogously, the terms quadratic in the top squarks are: 
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Figure 1: Diagrams that contribute to the two-loop effective potential and affect the 0{ata s ) 
calculation of the neutral Higgs boson masses. 

and the field-dependent entries of the stop mass matrix are, neglecting D-term contributions 
that vanish for g = g' = 0: 

m 2 L = m 2 Q + h 2 |F 2 °| 2 , m 2 R = m 2 u + h 2 \H 2 \ 2 , (13) 

X = \X\e^ = h t (A t H^ + ^H^ , (0 < <p < 2vr) , (14) 

where tuq, nifj and At are the soft supersymmetry-breaking mass parameters of the stop sector. 
We assume here \x and A t to be real, so that = (h%v%j\p2) \A t + /xcot/3| and {e tLp ) = 

sign {A t + //cot /3), but we do not make any assumption on their sign. 

The two phases tp and yp depend on the Higgs background. Therefore, in the evaluation of 
the derivatives of V e s, their contribution should not be neglected. To simplify the calculations, 
we choose to redefine the fields in such a way that the top and stop mass matrices become real: 

t' L = e-^t L , t' R = e^t R , t' L = e-^t L , t' R = e^t R . (15) 

This redefinition allows us to combine and t R , in the usual way, into a four-component Dirac 
spinor, with a real field-dependent mass mt = ht l-fffl- Moreover, the field-dependent stop mass 
matrix is now real and symmetric, thus it can be diagonalized by the orthogonal transformation: 




U)- 

The field-dependent stop masses are: 

(m 2 L + m 2 R ) ± \J (m 2 L - m 2 R ) 2 + 4|X| 2 



2 1 

nir = - 

*1,2 2 



(16) 



(17) 



and the mixing angle 0i is also a field-dependent quantity, defined by: 



2 \X\ 

sin2%= ' ' ■ (18) 
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Notice that, in this case, < 9^ < ir/2, in contrast with the usual field-independent definition 
for the angle 9^ that diagonalizes the stop mass matrix at the minimum, 



sin 29; 



2 rm (A t + fj, cot p) 



2 

ra- 
ti 



(19) 



which leads to —tt/2 < % < tt/2. The redefinition ([0]) has no effect on almost all field-dependent 
interaction vertices, with the only exception of the ones involving top, stop and gluino, which 
acquire a dependence on the phase difference (f — <p): 



-ttg 



V2g s (t L gTt L e^ 



tngTt R e 



+ h.c, 



(20) 



where T are the SU(3) generators in the fundamental representation and all color indices are 
understood. 

Before presenting the results for the 0{pna 3 ) corrections to the neutral Higgs boson masses, 
we discuss the general structure of the pure SQCD corrections to the one-loop 0(at) results, 
namely the 0(atd^) terms for generic n > 0. For the computation of this class of corrections, 
the effective potential can be expressed as a function of five field-dependent quantities, which 



can be chosen as follows. The masses m 2 



m~ and m~ , the last two as defined in 



i 



t 2 ' 



eq. ( |i~7| ) . The mixing parameter c 2 ^ = 1 — sin 2 29^ , where sin 29 ^ is given in eq. (|18|) . Finally, a 
parameter that, according to eq. (|20|), should be a function of the phase difference ip — tp: 
conveniently choose it as = cos (if — tp), where 



wc 



cos (ip — If) 



Re(X) Re(X) + Im(X) Im(X) 
\X\ \X\ 



(21) 



X and X are defined in eqs. (11) and (Q), respectively, and (c^) = ±1. A sixth parameter, the 
gluino mass rrig, does not depend on the Higgs background: we can restrict ourselves to positive 
values of rrig if we neglect possible CP-violating phases and we allow for arbitrary signs of [i and 
A t . 

According to eq. (||), in the limit of neglecting the SU(2)l x U(1)y gauge contributions 
beyond the tree level, the radiative corrections to the neutral CP-even Higgs boson mass matrix 
can be parametrized as: 

(M 2 3 ) ij = (Ml)°. + (AM 2 s )^, (22) 

where 
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(23) 



is fully determined by the input parameters mz, tua and tan/3 since, at 0(a s at), m 2 z = (g 2 + 
g' 2 )(v 2 + v 2 )/4:. The corrections that have not been reabsorbed in (M. 2 s ) a are contained in: 
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(24) 



Exploiting the field-dependence of the parameters, a wise although lengthy application of the 
chain rule for the derivatives of the effective potential leads to: 
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(AM* s ) u = ^ t V4^3, (25) 
(AM |) ^ = h 2 tf x m t s 2e F 2 + ^ h 2 t A t 11 4g (F 3 + AF 3 ) , (26) 
22 = 2 h 2 m 2 t F l + 2 ft 2 A t m t s 29 (F 2 + AF 2 ) + ^h 2 A 2 t s 2 e (F 3 + 2 AF 3 ) , (27) 



(AM 2 ) 



where s 2 g = sin2^ refers to the stop mixing angle defined in the usual field-independent way 
[see eq. (p^j)]. The functions Fi (i = 1, 2, 3) can be decomposed as Fi = F{ + AF{, where the AFi 
include the renormalization of the common factors multiplying F, (i.e. h 2 , mt, s 2 g), and 

~ _ d 2 V t d 2 V d 2 V 
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In the above equations, c 2S = 1 — s 2 g and the derivatives are evaluated at the minimum of V e g. As 
can be seen from eqs. (|25|)-(p7|), in every entry of AM. 2 S the Fi terms are multiplied by different 
combinations of \x and At- These two parameters do not renormalize in the same way, thus the 
contributions induced by their renormalization cannot be absorbed into the Fi, but should be 
separately taken into account. However, since at 0(a s ) [i does not renormalize, we have inserted 
in eqs. (26)-(|27|) only the two factors AF 2 and AF3 that take into account the renormalization 
of A t . 

To evaluate the functions (F\, F 2 , F 3 ), two strategies come to mind: i) evaluate explicitly the 
effective potential and then differentiate with respect to the relevant field-dependent quantities; 
ii) use a well-known fact, that the derivative of a bubble diagram with respect to an internal 
mass is still a diagram of the same kind, to compute directly the derivatives of the effective 
potential, without evaluating the effective potential itself. In our calculation of the 0(ata s ) 
corrections to the neutral Higgs boson masses we followed the latter strategy. The corresponding 
corrections to the effective potential can be found in eq. (4) of the second paper of ref. jf5|, with 
the understanding that the last two terms, proportional to mgmts 2 g, should be multiplied by 
c V {p. However, only the second derivatives of c^ with respect to the fields P, are different from 
zero at the minimum of the potential. Therefore, this extra term does not contribute to the 
expression for m^ in the decoupling limit of very large m&, where the results of that paper are 
applicable. 
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We give now the explicit expressions for the (D(ata s ) contribution to the functions (F\, F2, F3). 
For completeness, we recall first the one-loop result |H|: 
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where N c = 3 is a color factor. We assume that the 0{oct) one-loop contribution is written 
in terms of top and stop parameters evaluated in the DR scheme [v\ and t> 2 are automatically 
defined in the DR scheme by eqs. (|j) and (g), and the same is true for tan/3 = 1)2/ vi\. The two- 
loop 0{atQ. s ) contributions to the functions (F±, F2, F3), in units of g 2 Cf N c /(16 it 2 ) 2 (where 
Cf = 4/3), and in the DR renormalization scheme (here and hereafter DR quantities will be 
denoted by a hat), are given by: 



2( 



-6 1 



mi 



ln^| + 5 In 



2 2 

-M^+hi 



??? 



2 2 
2 h t 2 



mi 



/ m~ mf 

l n 2 _*i 1 i n 2 _J2 

v q 2 + Q 2 



s 29 



2 

mf 
n_ 

2 

*2 




In 



2 

mf 
2 ti 

m 



ti J 



+ /1 (mi , mj , m| , m| , s 2e , Q) + /1 (m t , m § , m| , m| , -s 2e , Q) , 



<2 



(32) 




m~ + m~ ra- 
ti *2_ 1 2 ti 

2 2 2 

mf — mf mf 

tl 12 t2 



+ S 2B 



mf 
tj_ 

2 

mf 

12 



1-ln 



Q 2 



mf 



mf 

t_2_ 

2 

*1 



+ / 2 K,^,mf i ,m|, S2e ,Q) 



5 In. 



m~ 



2 



' 2 , m | 2 , m | N 



m 



1-ln 



111 



to- 



rn 



f 2 {m t , m- g , m| , m\ ,-s 2e ,Q) , 



(33) 



16v^ 



F W (3 + 9c 2 e)+4 



3 + 13 c 2 , 



2 

mf 



tin 



ti 




+ 3 



2 2 
mf + mf 
*i 12 



m 



m 



ti 



hi 



2 

Q 2 



In 



2<' 

Q 2 



c 29 



/ 2 2 1 

mf + mf 

*1 *2 



m 



t-2- 



In 



2 

mf 
2 ti_ 

2 
mf 
£2 



2 2 
mf + mf 
ti 12 

(m~ — m? ) 2 

v *i 12 ' 



2 m i 
m?1 ln ~Q 2 



m 



ti 



"% 111 Q2 



111 



mf 

£l_ 

2 

mf 

12 



s 26» 



2 2 

mf m| 

£]_ 1 t_2_ 

2 ' 2 

mf mf 

12 *1 



7 



+2 In 



2 2 

mf mf 

*1 *2 



Q 4 

+ f 3 (mt,mg,m 



mf (m 



In 



*2 y 



<2 2 



+ 



4 

m~ 

52 



In 



9 

mf 

52 



mf 



'2 



m 



t2 



,Sz6,Q) + h{m t ,m~ g) m 



2 2 
fa'™*! 



Q 2 
S2e , Q) • 



(34) 



In the above expressions, Q indicates the DR renormalization scale, and the functions contain 
contributions coming from the top-stop-gluino diagrams (fig. |l]d): their explicit expressions are 
presented in appendix A. 

To obtain the 0{atOi s ) corrections to the Higgs mass entries, we also need explicit expressions 
for the AFi terms of eqs. (|26])-(p7|). In the DR scheme, and in units of g 2 Cf N c /(16 7r 2 ) 2 , they 
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A comment on eqs. (p2|)-(p6[) is in order. These equations show an explicit dependence on the 
renormalization scale Q, connected with our choice of expressing the top and stop parameters in 
the DR scheme. This dependence is cancelled by the implicit dependence of the DR parameters, 
so that the entries of AA^| are scale-independent. This fact becomes manifest if we reexpress 
the top and stop parameters in a physical scheme such as the on-shell (OS) scheme. To ensure 
this scale-independence, it was crucial to include the contributions induced by dV/dc^. If 
these terms were neglected, and the limit m^ — ► oo were taken, one would still find a scale- 
independent result for the (D(ata s ) corrections to m^ (thanks to the fact that, in such a limit, 
mh does not depend upon m^), but not for the corrections to m#. It may be useful to recall 
that, because v\ and vi are automatically defined in the DR scheme, (A^|)° has an implicit 
scale dependence, since Alp in eq. ( p^ ) contains tan/?. This residual scale dependence could 
be removed by including the momentum-dependent parts of the self-energies in the two-loop 
computation. 

To obtain the O(atot s ) corrections in some other renormalization scheme, R, we just need 
to shift the top and stop parameters appearing in the one-loop term. Indicating the general 
mass in the DR scheme as m DR , and the same quantity in the R scheme as m, we can write the 



one-loop relation m 
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8m. Then, once the one-loop contribution is written in terms of 



R quantities, the two-loop 0(at<y s ) corrections in the R scheme can be obtained through: 
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where all the quantities that appear in eqs. (37-41) are meant in the R scheme. The explicit 



,2 _2 

to the OS renormalization scheme, are listed in appendix B. 



expressions that allow us to perform the one-loop shift for (nit, m f 1 i m f 2 ' s 20) At), from the DR 



4 Discussion 

Using the formalism of the effective potential, we have obtained complete, explicit, analytical 
expressions for the 0(ata s ) two-loop corrections to the MSSM mass matrix for the CP-even 
Higgs bosons. Our input parameters are: (mz, tua, tan/3), already appearing in the tree-level 
result; the parameters of the top and stop sectors, appearing in the 0(a.t) one-loop correction, 
for example (nit, m~ , m| , s 2 #, a); the gluino mass and the strong coupling constant, appearing 
only at the two-loop level. We have presented our results in such a way that the input parameters 
of the top and stop sectors can be given either in the DR scheme or in some version of the on-shell 
scheme. Also, we have included in appendix C the 0(ata s ) corrections to the relation that gives 
m\ in terms of and tan j3. This result can be useful if one deals with models that predict the 
low-energy values of the soft supersymmetry-breaking parameters. 

Our effective potential calculation is equivalent to the evaluation of the Higgs self-energies 
in the limit of vanishing external momentum. A diagrammatic computation of the two-loop 
0(a,ta s ) contributions to the Higgs boson self-energies at zero external momentum has been 
performed in []T3| . Analytical formulae, valid in the simplified case of degenerate soft stop masses 
and zero mixing (with fi = At = 0), have been presented in the first paper of ref. |ll|. For 
arbitrary values of the top and stop parameters, however, the complete analytical result of 
is far too long to be explicitly presented, and is only available as a computer code [pl|. We have 
checked that, in the case of zero mixing and degenerate stop masses, our results coincide with 
those of [O]. Moreover, after taking into account the difference in the definitions of the on-shell 
renormalized angle 8^ (see appendix B), we find perfect agreement with the numerical results of 
1 18], for arbitrary values of all the input parameters. 

A calculation of both 0(ata s ) and O(o%) two-loop corrections to the lightest Higgs boson 
mass nih, based on the formalism of the effective potential, has been presented in Jl^]. In these 
papers, however, the dependence of the stop masses and mixing angles on the fields Pi (the 
CP-odd components of the neutral Higgs fields) is not taken into consideration. Therefore, 
the 0(ata s ) corrections to the input parameter niA are not evaluated. If one wants to relate 
the input parameters to measurable quantities, the computation is applicable only in the limit 
niA 3> niz, in which is nearly independent of niA- However, the results of [15] allow to 
express nih as a function of the (unphysical) renormalized parameter m| in the DR scheme. 
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Moreover, while an analytical formula for V e g is given (in terms of m\ and of the fields Si), the 
results of [15 1 for are available in numerical form, and simple analytical formulae are provided 
only in the case of universal soft stop masses and small stop mixing. 

The corrections controlled by a s introduce a new mass scale in the prediction of the MSSM 
neutral Higgs boson masses, namely the gluino mass m g . To avoid dealing with many different 
scales, m g is usually set to be of the same magnitude of the stop masses. Notice that, in 
a scenario where stops and gluinos are all heavy and approximately degenerate, with masses 
O(Ms), only the function F\ contains large logarithms of the ratio mt/Mg. Instead, as we 
have explicitly checked, F2 and F3 are finite in the limit mt — * 0. They contribute only to 
the matching conditions between the MSSM and the effective theory below the scale Ms, to 
be identified at 0(atct s ) with a two-Higgs-doublet version of the Standard Model. However, 
if the term At + //cot/? is very large, or and m,R are very different, or both, the two stop 
mass eigenstates can have a wide mass gap, and large logarithms of the ratio mi /mj can 
then be present in all the Fi terms. It should be mentioned that in this case the low-energy 
effective theory is different from a two-Higgs-doublet version of the SM, and indeed much more 
complicated already in the case of small stop mixing |J , not to mention the difficult case of large 
stop mixing. 

We can also envisage a scenario in which the gluino is much heavier than the top and the 
stops. Eqs. (|Al[) -([A^) contain terms proportional to powers of m g that can be potentially large. 
This powerlike behavior is actually cancelled in the OS schemes by the finite parts of the relevant 
shifts. However, as already noticed in [13], the gluino does not fully decouple, and increases 
logarithmically with m g when the latter becomes very large. On the other hand, it must be 
noticed that, in the DR scheme, some terms proportional to m g and m? are not cancelled, and 
in the limit of heavy gluino the two-loop corrections to the Higgs masses can become very large : 
this is related with the non-decoupling properties of mass-independent renormalization schemes 
such as DR. 

Prom eqs. (|A1|) - (|A3| ), we can derive approximate formulae for the two-loop corrections to M s 
in the case of large gluino mass, keeping the leading terms in an expansion of the complete result 
in powers of m g . Doing so, the asymptotic behavior is approached quite slowly as m g increases. 
The reason is that some terms that are formally suppressed by inverse powers of m g are indeed 
enhanced by large numerical factors. In order to get an accurate approximation to the correct 
result, it is then preferable to include also the next -to-leading terms in the expansion in powers 
of nig. Specializing for simplicity to the case of degenerate soft stop masses and fj, = A t = 0, so 
that = m~ 2 = m| and the only non-zero correction to the Higgs mass matrix is (A_A4|) 22 , 
we find: 

(AX!) = |4,n?f^-l-01n4-31„^ + 2>" 2:! il 
V • s /22 87r 4 1 y 3 mf mf m| J 

h 2 t9 2 s m 2 t [32tt 2 , n 2 2 . 160m?; + 112m? 352 2l m m? 

o4 7T 4 mi \ 3 3 mf mf 



224 m? +288 m? m? , _ „ x _ m? 



ln — § + 32 (m t 2 + m?) In 2 — f - 32 m\ In 2 — | J + O ( m^ 4 J . (42) 



m| m| 



m? \ / 
mf J V 



In deriving eq. (42), we have assumed an on-shell renormalization for the top and stop pa- 
rameters. As anticipated above, if we were to write the one-loop corrections in terms of DR 
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Figure 2: The mass mh from our effective potential calculation, as a function of the gluino mass 
rrig, for (a) no mixing (mj = mx = 1015 GeV, S20 = 0) and (b) large mixing (mj = 1175 GeV, 
mi 2 = 825 GeV, S20 = 1) in the stop sector. The other MSSM parameters are: tua = 500 GeV, 
tan/? = 10, /x = and m t = 175 GeV. The dashed line in fig. |2|a is the approximate result of 
eq. ©. 



parameters, we would find in AA4| terms proportional to m|. In fact, as can be seen from 
eq. (|B3|), the finite shift <5m? scales for large m g as m~ (ln(m?/Q 2 ) — 1), and cancels in the OS 
scheme similar terms present in the fi. 

Some representative results for the heavy-gluino limit are shown in fig. ||[ We have plotted 
mh, as obtained from our complete formulae, as a function of the gluino mass. We employ on- 
shell top and stop parameters (for their definitions, see appendix B). We consider two examples, 
with degenerate soft stop masses and either no mixing or large mixing in the stop sector. The 
numerical inputs we use are: tua = 500 GeV, tan/3 = 10, fj, = 0, mt = 175 GeV, and the values 
raj = mi 2 = 1015 GeV, S20 = (no mixing, fig. ^a) and = 1175 GeV, m^ 2 = 825 GeV, 
s 26 = 1 (large mixing, fig. Qb). For the one-loop 0(ctt) corrections to M. 2 S we have used the 
effective potential result of ||. As can be seen from figure m^ reaches a maximum at low values 
of mg, decreases for intermediate values and then increases logarithmically when mg becomes 
very large. Comparing the cases of no mixing and large stop mixing, we see that in the latter case 
the peak of m^ at small mg is much more pronounced, and the asymptotic increase of m^ starts 
at higher values of mg. In fig. |2|a we have also plotted m^ as obtained with the approximate 



formula of eq. (42), valid in the case of no mixing. One can see that eq. ( [42]) approximates very 
well the complete result when mg > 2 m^ . 
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Appendix A: Analytical expressions for the gluino contributions 

The explicit expressions of the contributions to the functions F 2i (i = 1,2,3) coming from the 
top-stop-gluino diagrams (fig. [j]d) are rather long but, as apparent from eqs. (|32"|)-([34|), they 
possess useful symmetry properties under the exchanges m~ «-> mf , s 2 q <-> — s 2 q. In units of 
g 2 Cf N c /(16tt 2 ) 2 , where Cf = 4/3 and N c = 3 are color factors, the explicit expressions of the 
functions fi (i = 1,2,3) are: 
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where A = m\ + mf + mf — 2 (m? m? + 



,2 ^,2 



and the function <I> is defined as in 



2(1: 
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$ (x, y, z) = — ^2 lnx + lnx_ — Inu Inv — 2 (Lia x + + Li2 x_) + — j , (A4) 

where the auxiliary (complex) variables are: 

u = , v = -, A = ^(1-u-v) 2 -4uv , x± = - (1 ± (u - v) - A) . (A5) 

The definition ( |A4| ) is valid for the case x/z < 1 and y/z < 1. The other branches of $ can be 
obtained using the symmetry properties: 

$(x,y,z) = $(y,x,z), x$(x,y,z) = z$(z,y,x). (A6) 



Appendix B: Shifts of the parameters to the on— shell scheme 

In the OS renormalization scheme, the masses of all particles are defined as the poles of the 
corresponding propagators. As an example, for a scalar particle with squared mass ra 2 the 
relation between the DR and OS definitions of the mass is: 



bra 



(m 2 ) DR - K) 



2nOS 



Ren(m 2 ) 



(Bl) 



where li(m 2 ) is the finite part of the self-energy of the particle, evaluated at an external mo- 
mentum equal to the mass itself. In the following we list the shifts to the on-shell scheme for 
the top and stop masses []: 
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ti 



bra 2 
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[(1,820) <- (2,-52*)] , 



(B3) 



(B4) 



where the notation (1, S20) (2, — S20) i n eq.(B2) means a term that is obtained from the previous 
ones inside the square bracket with the exchange raf <-> m| and the replacement S20 — ► —s 2 e- 



ti ""ta 

Formulae for the DR— OS shifts, specialized to the case 



1 Similar results have been presented in pl[ jnj. 
= rrlji, can also be found in However, we disagree with JTsJ on the contributions to the stop self-energies 
coming from the diagrams that involve the 0(a s ) four-squark vertex. 
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The notation of Eq. (B4) implies that /m| can be obtained from the right hand side of 



eq. QB3D with the above substitutions. The quantity 5 that appears in eq. QB2| ) is a constant 
that depends on the regularization. In dimensional regularization 5 = —4, while in dimensional 
reduction 5 = —5. In eqs. (|B2|)-(p4]), Bq denotes the finite part of the Passarino-Veltman 
function, i.e. : 



dx In 



'1 — x) m\ + x m\ — x (1 



x ) p 



IE 



Q 2 



(B5) 



An explicit expression for Bq can be found e.g. in [22|. 

Due to the relative freedom in the choice of the renormalization conditions for the stop sector, 
several definitions are possible for the shift in the mixing angle 9^ (for a discussion on this point, 
see [23 1 and references therein). We choose the following "symmetrical" definition: 
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where lii 2 (p 2 ) is the off-diagonal self-energy of the stops: 
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(B7) 



Finally, taking into account that fi and tan/3 do not get any O (a s ) correction, the shift for the 
soft term A t can be easily derived from the (field-independent) definition of S20, eq. (|l9|): 



5 At 



ti 



5m i 2 Ss2e_ _ 5m t 



m~ 



S29 



(A t + n cot /3), 



(B8) 



where 5s2e/s2e = 2 cot 29^59^. Eq. ( |B8[) can be treated as a definition of At in our on-shell 
scheme. 



Appendix C: Two-loop corrections to the CP-odd Higgs mass 

We present here the two-loop 0(oeta s ) corrections to the mass of the CP-odd Higgs boson, A. As 
apparent from eq. (^), m\ depends on the value of the soft super symmetry-breaking parameter 
7713, thus we included it among the MSSM input parameters. However, a calculation of the loop 
corrections to the relation between m\, m| and tan/3 may be useful for discussing models that 
predict the values of the soft super symmetry-breaking masses, and in particular m 2 . 

Starting from eqs. (|6|) and (^), and following the same line of reasoning as in Section 3, we 
find that: 



^-iE^l^ + jra-^J- (C1) 

where niA is the Q{a.ta™) loop-corrected A mass evaluated in the effective potential approach, 
i.e. neglecting the corrections that depend on the external momenta. The function Fa can be 
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decomposed as Fa = Fa + AFa , where AFa contains terms coming from the renormalization 
of the parameters that multiply Fa in eq. (|Cl|), and Fa is defined as: 

dV 2\i cot(3{s 2 e)~ 2 dV 



F A 



dV 



dV 



dm 2 dm 2 

12 tl 



+ 



m i - m i dc 2§ 



A t {m 2 h -m 2 i2 ) dcp$ 



(C2) 



The definitions of the field-dependent parameters m~ , m 2 ^ , c 2 ^ and c^q, are given in Section 
3. The one-loop 0{at) contribution to Fa is known 
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Assuming a DR renormalization for the parameters ht , ^ , m 2 and m-^ that enter in eqs. ( |C1| ) 



and (C3), the two-loop 0{ata s ) contribution to Fa in the DR renormalization scheme is given, 
in units of g 2 C F N c /(16ir 2 ) 2 , by: 



hit 



16 vr 2 



F 



8 - s 2 e 2 



I'll; 



mf + m 2 m? 

tl *2_ 1 ti_ 

2 2 2 

mf — mf mf 

tl t2 12 ■ 



9 

mf 



+2 ml In 2 -A - m? In 2 + 



*2 



Q 2 



W; — m, 

11 12 



/ m| m| \ 2 



+ fA(m t ,m § ,ml,m 2 2 ,s 2e ,Q) - f A (m t , m s , m\ , m\ , -s 26 , Q) 



(C4) 



The function Ja contains contributions coming from the top-stop-gluino diagrams (fig. [yd), and 
its explicit expression in units of g 2 Cp N C /(1Q 7r 2 ) 2 is: 
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where A and y, z) are defined at the end of appendix A. 
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